A banded preconditioner for the two-sided, nonlinear space-fractional diffusion equation
نویسندگان
چکیده
The method of lines is a standard method for advancing the solution of partial differential equations (PDEs) in time. In one sense, the method applies equally well to space-fractional PDEs as it does to integer-order PDEs. However, there is a significant challenge when solving space-fractional PDEs in this way, owing to the non-local nature of the fractional derivatives. Each equation in the resulting semi-discrete system involves contributions from every spatial node in the domain. This has important consequences for the efficiency of the numerical solver, especially when the system is large. First, the Jacobian matrix of the system is dense, and hence methods that avoid the need to form and factorise this matrix are preferred. Second, since the cost of evaluating the discrete equations is high, it is essential to minimise the number of evaluations required to advance the solution in time. In this paper, we show how an effective preconditioner is essential for improving the efficiency of the method of lines for solving a quite general two-sided, nonlinear space-fractional diffusion equation. A key contribution is to show how to construct suitable banded approximations to the system Jacobian for preconditioning purposes that permit high orders and large stepsizes to be used in the temporal integration, without requiring dense matrices to be formed. The results of numerical experiments are presented that demonstrate the effectiveness of this approach.
منابع مشابه
Preconditioned Generalized Minimal Residual Method for Solving Fractional Advection-Diffusion Equation
Introduction Fractional differential equations (FDEs) have attracted much attention and have been widely used in the fields of finance, physics, image processing, and biology, etc. It is not always possible to find an analytical solution for such equations. The approximate solution or numerical scheme may be a good approach, particularly, the schemes in numerical linear algebra for solving ...
متن کاملThe new implicit finite difference scheme for two-sided space-time fractional partial differential equation
Fractional order partial differential equations are generalizations of classical partial differential equations. Increasingly, these models are used in applications such as fluid flow, finance and others. In this paper we examine some practical numerical methods to solve a class of initial- boundary value fractional partial differential equations with variable coefficients on a finite domain. S...
متن کاملAn Implicit Difference-ADI Method for the Two-dimensional Space-time Fractional Diffusion Equation
Fractional order diffusion equations are generalizations of classical diffusion equations which are used to model in physics, finance, engineering, etc. In this paper we present an implicit difference approximation by using the alternating directions implicit (ADI) approach to solve the two-dimensional space-time fractional diffusion equation (2DSTFDE) on a finite domain. Consistency, unconditi...
متن کاملOn the split-step method for the solution of nonlinear Schr"{o}dinger equation with the Riesz space fractional derivative
The aim of this paper is to extend the split-step idea for the solution of fractional partial differential equations. We consider the multidimensional nonlinear Schr"{o}dinger equation with the Riesz space fractional derivative and propose an efficient numerical algorithm to obtain it's approximate solutions. To this end, we first discretize the Riesz fractional derivative then apply the Crank-...
متن کاملFourth-order numerical solution of a fractional PDE with the nonlinear source term in the electroanalytical chemistry
The aim of this paper is to study the high order difference scheme for the solution of a fractional partial differential equation (PDE) in the electroanalytical chemistry. The space fractional derivative is described in the Riemann-Liouville sense. In the proposed scheme we discretize the space derivative with a fourth-order compact scheme and use the Grunwald- Letnikov discretization of the Ri...
متن کاملذخیره در منابع من
با ذخیره ی این منبع در منابع من، دسترسی به آن را برای استفاده های بعدی آسان تر کنید
برای دانلود متن کامل این مقاله و بیش از 32 میلیون مقاله دیگر ابتدا ثبت نام کنید
ثبت ناماگر عضو سایت هستید لطفا وارد حساب کاربری خود شوید
ورودعنوان ژورنال:
- Computers & Mathematics with Applications
دوره 66 شماره
صفحات -
تاریخ انتشار 2013